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Abstract 

We compute the four-loop QCD contribution to the electroweak p parameter 
induced by the singlet diagrams of the Z-boson self-energy. The numerical impact 
on the weak mixing angle and the VK-boson mass is small. 

PACS numbers: 12.38.-t, 14.65.Ha, 13.66.Jn 



1 Introduction 



The electroweak p parameter as introduced by Veltman pQ measures the relative strength 
of the charged and neutral current. Considering QCD corrections it can be written as 

P = 1 + Sp, (1) 

with 

n zz (o) iw(o) 
6p = ~Ml W (2) 

rizz(O) and II^v^(O) are the transverse parts of the W- and Z-boson self-energies evalu- 
ated for vanishing external momentum. The parameter bp enters a variety of quantities 
which are determined from experiment with an enormous precision. In particular, it en- 
ters the relation between the iy-boson mass, My/, the fine structure constant, a, the 
Fermi constant, Gp, and the Z-boson mass, M z , which is given by [2] 



Ml- = ^11 



Ana 



\\ V2M Z G F {1 - Ar) 
The quantity Ar is conveniently parameterized in the form 



(3) 



Ar = Aa - ^fbp + Ar rem , (4) 
s w 

with cw = M w /M z and = 1 — c^. Aa contains contributions from light fermions 
giving rise to a correction of about 6%. The leading corrections proportional to GpM^ are 
incorporated in bp and amount at one-loop order to roughly —3% whereas the remaining 
part is small. 

Eqs. (JHJ) and (pEJ) can be used to predict My/ where the formula 



SM W = ^^^rbp, (5) 
c w s w 



immediately accounts for the dominant shift in My/ due to the corrections to the p 
parameter. We can also look at the change of the effective leptonic weak mixing angle, 
sin 2 ^^*, defined through the coupling of the Z-boson to leptons. The leading universal 
corrections originating from bp can in analogy to Eq. (JHJ) be written as 



5sin 2 ^ ff pt = - f wS \ bp. (6) 
z w s w 



Currently the uncertainties for My/ and sin 2 O 1 ^ are given by bMyy = 34 MeV and 
^sin 2 ^^* = 1.7 x 10 -4 J3], respectively. However, a future linear collider running at the 
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Figure V. Sample three- and four- loop singlet diagrams contributing to the p parameter. 
In the fermion loops either top- or bottom-quarks are present. 

Z-boson pole, the so-called Giga-Z option, and around the W-pair threshold might reduce 
the uncertainties to 5M W = 6 MeV and 5 sin 2 #' e ff pt = 1.3 x 10" 5 @]. 

The one-loop corrections to p have been computed in 1977 pQ and also the two-loop 
QCD corrections are known since almost 20 years Roughly 10 years ago the order 

GVM 2 a 2 QCD corrections [HUH] constituted one of the first applications of the three-loop 
massive vacuum integrals. At three-loop order for the first time a new kind of Feynman 
graphs has to be considered, the so-called singlet diagrams as shown in Fig. [T] which only 
contribute to the Z-boson self-energy. They are characterized by the fact that in contrast 
to the non-singlet contribution the external Z-bosons couple to different fermion lines. We 
want to note that the singlet contribution forms a finite and gauge independent subset. 
At three-loop order it completely dominates the numerical corrections if the MS definition 
is adopted for the top-quark mass. In the case of the pole mass definition the singlet part 
still amounts to about 30% of the total three-loop contribution. We want to mention that 
also two-loop [TUJ HI] and three-loop mixed electroweak/QCD [T2] and even three-loop 
pure electroweak corrections ^2] have been evaluated. Recently also corrections in the 
large Higgs boson mass limit have been considered [T3"lll4j. For non-universal corrections 
to Mw and sin 2 9^ we refer to [T5| ITflj. 

In this letter we consider the four-loop contribution to the p parameter originating 
from the singlet diagrams. In Fig. some sample diagrams are shown. This constitutes 
one of the first applications of the four-loop vacuum master integrals evaluated recently 
in Ref. [IZj. 

2 Technicalities 

Since the boson self-energies have to be evaluated for zero external momentum and only 
QCD corrections are considered, only the axial-vector part of the Z-boson correlator 
gives a non-zero contribution. Whereas for the non-singlet contribution the naive anti- 
commuting definition of 75 can be adopted, special care has to be taken in the singlet case. 
Actually, the definition of 't Hooft and Veltman JH] has to be adopted and additional 
counterterms have to be introduced in order to ensure the validity of the Ward identities. 
In the practical calculation we follow Ref. [1.9] and perform the following replacement in 
the axial- vector current 

tV = \e^° lvlplu . (7) 
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We pull out the e-tensor from the actual integral and consider instead the completely 
antisymmetrized product of the three 7-matrices which can be written as 

7 [1 W ] = ^( 7 Y7 ff -7 ff 7Y). (8) 

As a consequence we have to deal with an object with six indices. Thus, for zero external 
momentum we obtain 



Wzz - —r^zz 



H[fpcr] [v'p'tj 1 



144 

Ijt^H (Q) 

24 1 >H • 1 j 

In the practical calculation we consider the object II ^ for which we also perform the 
renormalization as described in the following. Thus, in Eq. © the limit D — > 4 has been 
considered where D = 4 — 2e is the space-time dimension. 

The additional finite counterterm is only needed to one-loop order, since the singlet 
diagrams appear the first time at three-loop level. For each axial- vector vertex a factor [201 

m 

Z« = 1 - C F ^ + 0(a 2 s ) , (10) 

7T 

with Cf = (N% — 1)/(2N C ) has to be considered. Furthermore, we have to consider the 
one-loop counterterms for the strong coupling constant and the top-quark mass defined 
by 

a° s = Z aa a s , m° = Z m m t , (11) 

where m t = mt(n) is renormalized in the MS scheme. The renormalization constants are 
given by 

Z as = l + \(-T 2 C A + \Tn f )^ + 0{al), 

Z m = l-^- + 0(« s 2 ), (12) 

with Ca = N c and T = 1/2. rif = 6 is the number of active flavours. The transition to 
the pole mass is achieved via 



mt(jj) 



/i 2 \ a s 



- 1 -4 1 "w 



M, . (13) 



We generate the Feynman diagrams with QGRAF [21 j and adopt with the help of the 
packages q2e and exp ESj the topologies and notation to the program performing the 
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reduction of the four- loop vacuum diagrams jz4j. As an output we obtain the corrections 
to the p parameter as a linear combination of several master integrals. All of them have 
been computed in Ref. [TTj . 

It is interesting to note that some of the master integrals are multiplied by spurious 
poles of order 1/e 2 . As a consequence, for these the (9(e) and even the C(e 2 ) contribution 
is needed. In the case of the master integral BB4 (which is the four-loop sunset vacuum 
bubble with one massless and four massive lines, see Eqs. (4.8) and (6.36) of Ref. [T7j) it 
happens that the coefficient of order e, which originally has only been evaluated in numer- 
ical form enters the pole part of 5 p. Thus an analytical expression can be deduced 
(in Eq. (fT4*|) below denoted by BBA^) which perfectly agrees with the known numerical 
result. Furthermore, we have obtained an analytical expression for the coefficient of order 
e 2 , by combining the numerically known value with the basis of transcendentals known 
from an independent investigation |25j . These two coefficients read 



(n-4) 



BB4 = J 4 ^e n BB4 

n>0 

m 1976975 1792 >/oN 

BB4 (1) = + C(3) 

7776 9 SV ; 

= -14.897726533029588869214274870082319534267 . . . , 

(2) 72443143 47488 >/oN 8704 

BB4 (2) = + C 3) C(4) 

46656 27 SV ; 3 SV ; 

1024, 4 2048 > , 2 8192 
+ — In 4 2 - —C(2) In 2 2 + — a 4 

= -1678.886929107772963403030310267917509151 . . . , (14) 

where J is the one- loop tadpole, ((n) is Riemann's zeta function and 

a 4 = Li 4 (l/2) « 0.51747906167389938633. (15) 

Let us mention that we performed the calculation using an arbitrary gauge parameter 
of the QCD gluon propagator, £. As expected the final result is independent of £ even 
before inserting the values for the master integrals. This constitutes a nice check of our 
result. 



3 Results and discussion 

Let us in the following present our analytical result and discuss its numerical implications. 
For completeness we also repeat the QCD corrections up to three-loop order. For the MS 
definition of the top-quark mass we obtain 



Sp 



MS 



Air 



+ 



4tJ 



26459 25064 



51 



C(2) 



3560 
~27~ 



C(3) 



1144 , s 
+ — C(4) 



16 



B A - -D 3 + 882S 2 + n f 
9 



50 112 „ s 64 w ^ 



5 



9 9 



mt 



+ 



An) 



^-4528C(3) 



20816 



2624 



5248 



In 4 2 + C(2)ln 2 2 



20992 



a 4 - 1232C(3)lnAr 

3 mi 



C(4) 



+ . . . 



(16) 



where the "—56^(3)" in the third line stems from the three-loop singlet diagram. In order 
a 3 only the singlet contribution is presented. Furthermore, we have 



G F m\ 
8^/2' 

^=Im(Li 2 (e m/3 )) w 0.26043413763216209896 

2 1 ^ 

16a 4 + -ln 4 2-4C(2) hi 2 2- y C(4) 
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6C(3)-^C(4)-6[Im(Li 2 (e^ 3 ))] 2 



-1.7628000870737708641 
-3.0270094939876520198 . 



In Appendix [21 we present the three- and four-loop result for the singlet contribution 
corresponding to Eq. (fT^j) retaining, however, the colour factors Cp, Ca and T. With the 
help of Eq. ()13j) one obtains the singlet result in the on-shell scheme. Together with the 
non-singlet terms one gets 



Sp 



OS 



3X t 1 + 



47T 



8 16 



C(2) 



+ 



16 



a. 



-V 



Ait J 



314 26504 



C(2) 



3416 



64 1144 
- — C(2) In 2 + -g-C(4) - - -D 3 + 882S 2 + n f 



27 
8 208 
9 



C(3) 
C(2) 



9 



64 



38 176 



-C(3) ) - 56C(3) + ( -- - — C(2) +n f (- + -((2) ) ) In ^ 



16 32 



Ml 



+ 



47T 

2624 
"9" 



256 11792 20816^ 



9 3 
In 4 2 + —((2) hi 2 2 



3 

20992 



., a 4 - 784C(3)ln^ 



+ . . . 



(17) 



with X t = G F M?/(8rr 2 V2). 

Inserting the numerical values for the constants in Eqs. ()16)) and (|T7|) and adopting 
p = m t and p = M t , respectively, the numerical corrections read 



MS 



bp 

5p os 



3xt 
3X t 



ry / qj 

0.19325— + (-4.2072 + 0.23764) — 

V 7T 



7T 



1 - 2.8599— + (-4.2072 - 10 

7T 



2 /a ^ 3 

3.2866 [ — 

7T 



+ 7.9326 



ex. 



71 
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where the three-loop contribution is split into the singlet (first number in round brackets) 
and the non-singlet piece. If we furthermore adopt a s (m t ) = 0.108 and a s (M t ) = 0.107, 
the expression for 5p looks like 

S p m = 3x t (l - 0.00664 - 0.00469 — 0.00013^ , 

5p os = 3X t U - 0.09741 - 0.01693 + 0.00031 J , (19) 

where the n th term inside the round brackets corresponds to the contribution of order 
GpM'^a'f 1 ^. One observes that the new four-loop singlet contribution is numerically 
small and amounts to about 3% of the three-loop result in the MS scheme and to less 
than 2% for on-shell top-quark masses. Note that the correction is positive in the on-shell 
and negative in the MS scheme. In the on-shell scheme the shift in M w and sin 2 9^ 
according to Eqs. (0) and © amounts to 0.175 MeV and 10~ 6 , respectively, which is 
significantly below the recent estimates of higher order contributions and variations of 
input parameters (13 HE] . 

It is interesting to mention that at three-loop order the singlet contribution completely 
dominates for MS top-quark masses and amounts to about 30% in the on-shell scheme. 
Thus, in case the same pattern also holds at four-loop order, the complete QCD corrections 
would be well under control. However, the numerical values in Eq. (fTHj) suggest that for 
some reason the four-loop singlet contribution seems to be accidentally small. 

Let us also comment on the dependence of the singlet contribution on the renormaliza- 
tion scale p which can be done separately from the non-singlet part. The latter is discussed 
in Ref. [§| (cf. Fig. 1 of Ref. [Hj). As far as the singlet contribution is concerned one ob- 
tains for the quantity (<5p os /(3X<) _ l) singlet the values {-0.00457, -0.00437, -0.00455} 
corresponding to p = {M t , M t /2, 2M t }. The /i-dependence, being formally of higher or- 
der, is less then 5% of the sum of the three- and four-loop singlet part which can be used 
as an estimate of the 0{af) term. 

In the remaining part of this section we briefly compare the numerical effect of the 
new terms with known corrections to 5p. In the on-shell scheme the three-loop QCD 
corrections of order a 2 s X t lead to a shift of about —10 MeV in the IV-boson mass and to 
+5 x 10~ 5 in the effective weak mixing angle. For Higgs-boson masses between 200 GeV 
and 300 GeV the three-loop corrections of order a s X 2 have the opposite sign and with 
roughly half the magnitude they are still relevant for the precision to be reached at the 
Giga-Z option of a future e + e~ linear collider. However, the pure electroweak corrections 
of order Xf are very small and give rise to corrections well below 1 MeV for the shift in the 
PF-boson mass. The same is true for the four-loop QCD singlet contributions considered 
in this letter. 

In conclusion, we computed the four-loop singlet contribution to the p parameter 
which constitutes one of the first applications of the four-loop massive vacuum integrals 
to a physical quantity. The numerical size of the corrections turn out to be surprisingly 
small and lead to a shift in the jy-boson mass below 1 MeV and to the effective weak 
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mixing angle below 10 -5 — beyond the accuracy forseen in a future Linear Collider. This 
illustrates the good convergence properties of the perturbation theory and confirms the 
stable predictions based on the three-loop corrections. However, for a definite conclusion 
also the non-singlet contribution has to be evaluated. 
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A Singlet contribution to the p parameter 

In this Appendix we present the three- and four-loop singlet result expressed in terms of 
C A = N c , C F = (A" 2 - 1)/(2JV C ) and T = 1/2. Furthermore, we keep the label n t which 
counts the number of massless quarks. The three-loop term can also be found in Ref. 



2 



= 3,-, ( C F T\ - 84C(3) + g 



C F \ - 336C(3) + 2400C(4)-1281n 4 2 



,2 



7064 



+ 768C(2) In 2 2 - 3072a 4 - 504C(3) In !-= + C A [ C(3) 

mf J V 3 

+ 3056C(4) - — In 4 2 + 640C(2) In 2 2 - 2560a 4 - 616C(3) In At 
3 mf 

+ n{T (™C(3) - 784C(4) + y In 4 2 - 128((2) In 2 2 + 512a 4 



+ 224C(3)ln g) +r (f-lf2 f( 3 )+224C( 3),f 



(20) 
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